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We study the spectral properties of infinite rectangular quantum graphs in the presence of a 
magnetic field. We study how these properties are affected when three-dimensionality is considered, 
in particular, the chaological properties. We then establish the quantization of the Hall transverse 
conductivity for these systems. This quantization is obtained by relating the transverse conductivity 
to topological invariants. The different integer values of the Hall conductivity are explicitly computed 
for an anisotropic diffusion system which leads to fractal phase diagrams. 



I. INTRODUCTION 



Quantum graphs have been the focus of much interest during the last thirty years [H El E]. These models which 
describe the propagation of a quantum wave within an arbitrary complex object are extremely versatile allowing the 
study of various interesting quantum phenomena. Quantum graphs appear in various fields such as solid state physics, 
quantum chemistry, chaology and wave physics. Basically, quantum graphs describe wave propagation through fine 
structures. In the field of quantum chemistry, graphs have been used to represent 7r-electronic orbitals in organic 
molecules formed with double chemical bonds In nanotechnology as for future quantum computer devices, they 
modelize thin conductor circuits that propagate information. They also describe fine superconducting circuits and 
wave guides leading to acoustic, optic and electromagnetic applications. Generally speaking, graphs constitute useful 
models for the description of quantum transport on connected systems W. 

In the context of quantum chaology, graphs have been the vehicle to confirm important conjectures about chaos 
signatures. Kottos and Smilansky discovered that periodic-orbit theory exactly applies to quantum graphs and that 
their spectra may obey the Wigner level spacing statistics under certain conditions [5]. A minimum of three incom- 
mensurate bond lengths is required for Wigner level repulsion to already manifest itself Furthermore, graphs are 
simple models of quantum scattering and a semiclassical bound has been obtained on their quantum lifetimes [7] . In 
the classical limit, the scattering is typically stochastic at each vertex of the graph so that the classical dynamics may 
be chaotic with a positive Kolmogorov-Sinai entropy per unit time on graphs with more than two vertices [8]. 

In this paper, we propose to study the integer quantum Hall effect (IQHE) on the basis of quantum graphs. This 
phenomenon which is manifested by the quantization of the Hall conductance is the object of many works since its 
discovery by von Klitzing [S]. Thouless et al. [TO] showed the important link between the Hall conductance and the 
energy spectra of independent electron's model, a result which led Osadchy and Avron to draw the phase diagram 
for Hofstadter's model |T5]. This phase diagram is fractal and depicts infinitely many phases, each one characterized 
by the integer value of the Hall conductance. Such fractal phase diagrams have also been studied in cold atomic 
systems submitted to artificial gauge fields [13] and in models featuring continuous potentials which exhibit classical 
chaos [mill]. 

In this context, quantum graphs are interesting because the propagative medium is formed of one-dimensional 
continuous bonds, in contrast to the fully discrete lattices of the usual tight-binding models. Yet, quantum graphs are 
simple enough that they can be investigated by analytical developments, as we shall show. 

We first study the energy spectra for infinite rectangular graphs perturbed by a strong magnetic field. We obtain 
the spectral properties of two-dimensional {ID) and three-dimensional (31?) systems. This allows us to compare these 
systems from the viewpoint of chaology. We show that the eigenvalue equation of quantum graphs can be mapped onto 
the generalized Harper equation in the case of 2D rectangular lattices. We then compute the transverse conductivity 
of the 2D system with Kubo's formula which describes the linear response of this system to an external electric field, 
and obtain the quantization law through topological arguments. We then show the link between the Hall conductivity 
and the energy spectra by computing this quantity for Fermi energies located inside the many gaps of the spectra. 
We conclude this work with the presentation of a fractal phase diagram which describes the IQHE on 2D as well as 
3D quantum graphs. 

The paper is organized as follows. The spectral properties of 2D and 2>D quantum graphs are respectively presented 
in Sees. [H] and HI The Hall conductivity on quantum graphs is studied in Sec. IV The conclusions are drawn in Sec. 
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II. TWO-DIMENSIONAL QUANTUM GRAPHS 



A. The eigenvalue equation 



A quantum graph is a set of vertices connected by unidimensional bonds on which a quantum wave propagates. 
Each bond of a rectangular lattice is characterized by a vertex coordinate (n, m) and by a direction b = 1,2, and will 
be labeled (n,m,b). Schrodinger's equation is satisfied on each bond 



dn, 



(1) 



where A'^"'™'''-' is the potentiel vector component along the bond b, k ^ \/2ME/h is the wave number, M and E 
are respectively the mass and the energy of the particle. In the following, we use units where M — 1 and h = 1, 
except otherwise stated. We consider an infinite graph forming a rectangular lattice submitted to a magnetic field 
B = 'V X A and we work in the Lorentz gauge A = {0,Bx,0), which maintains a discrete translational invariance 
along the y direction. We suppose that the different bonds are directed along the x (resp. y) axis with constant length 
Ix (resp. ly). The wave function ip'-"-'"^'^\x) (resp. ^'■"■'"^'^\y)) is defined between vertex {n,m) and (n + 1,to) [resp. 
(n, m) and (n, m + 1)]. Accordingly, we have that x € [0; Ix] and y £ [0; ly]. 
The solutions of Schrodinger's equation ([l]) are written as 



(2) 



We consider an anisotropic scattering model which results from the modification of the usual boundary conditions 
imposed at each vertex. We introduce two parameters s and r in these conditions 

^(«,m,i)(2. = 0) = = Ix) 

(3) 
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(4) 



where the particular case r — s = I corresponds to an isotropic graph. One has to impose the unitarity of the scattering 
matrix S. This matrix which plays an important role in the theory of quantum graphs links the in-components x'" 
to the out-components at each vertex. These components are written as follows for our model 



/ ^(n,m,l) 

^(n,m.2) 
'-■1 

^{n-l,m,l) ^^iki^ 



) 



\ 



^(n,m,2) 
^2 



(5) 



where (j')^^'^^ are the coefficients of the functions ^/;("'™''') as defined in Eqs. ([2]). The scattering matrix is then given 

by 



5* = 




(6) 
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and the unitarity condition S^S = SS^ = I sets r = s. The unitarity condition guarantees the conservation of 
probabihty at each vertex. This latter is expressed as 



j(n-l,m,l) _|_ j{n,m-l,2) _ j(n,m,l) _|_ j{n,m.,2) 

with the probabihty current densities : 



J 



{n,m.b) 



2i 



(n,m,6)* 



n,m,fc) 



J„l.{n,m,b)* \ 

dn J 



One can then write the solutions ^ in the following form 

,l,{n,m,i) ( \ _ f sin k{lx - x) sin kx 

sm klx sm klx 

r sm r sm kly 

where the coefficients fn,m are defined by 



m,2) ^(n,m,2) 
2 



(7) 



(8) 



(9) 



(10) 



Here, $ = Bl^ly/{2TT) is the magnetic flux through a unit cell and the condition ^ has been applied. 

The system being invariant under discrete translations along the y axis, we are led to consider new wave functions 
y(n,f>) defined in terms of the solutions ^("''"''') of Eq. ([!]) according to 



Similarly, the coefficients gn are defined as 

9n — fn^m ^ 

According to Eqs. (|9|, the new wave functions are written in terms of new coefficients (13 1 as 



(11) 
(12) 

(13) 



(ni)f ^ sinfc(;2,-a;) sinfca: ; 



sin kL 



9n+l 



sin klx 



The probability conservation Q implies 



sin fc^y sin fc/j^ 



2 sin fcL 



+ e ffn-i + ^ . , ,^ cos(27r$n - fcj^Zy) g„ 

sm Ktj, 

cos kL. sin fcZ 



= 2 cos fc/ 



If we set 



f ^ sin fc/, 



1 sin klx 



9n 



cos kly sin kl^ 

£ = 2 \ cos fcL ' 



sin fcL 



one finds the generalized Harper equation 



6'''"='"= gn+i + e ^'^-'^ Sn-i + 2Acos(27r$n - /c^^^j,) gn = £ gn 



(14) 
(15) 



(16) 

(17) 
(18) 

(19) 
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which reduces to Harper's equation for Hofstadter's model when = ly and r = 1. This shows that the problem of 
the quantum graph can be mapped onto a similar problem for the anisotropic tight-binding Hamiltonian with the 
transfer coefficients ta and tb [T7]. The correspondence is estabhshed with the anisotropy ratio A = th/ta- However, 
an important difference is that the parameter A and the energy £ are independent in the tight-binding model, albeit 
they both depend on the wave number k — ^/2E in the quantum graph. In this regard, the eigenvalue problem is 
more complicated in quantum graphs. 

We notice that we would have obtained the dual generalized Harper equation 



" h. 



2 E 

m+i + e~''^y^^ hm-1 + J cos(27r$m + kj^) /i„ = - /i„ 



(20) 



if we had used another gauge with the vector potential A' = {—B y, 0, 0). In principle, the eigenvalue equations p9| and 
(20 1 should lead to the same energy spectrum since they correspond to the same magnetic field B = 'VxA' = 'VxA. 



Indeed, Eq. (20 1 is derived from Eq. (19) thanks to the duality transformation 



+ 00 

E 



gn e 



(21) 



Our aim is to obtain the spectra in the plane of the magnetic flux <f> versus the wave number k — \/2E. This latter 
is directly related to the energy E = fc^/2 and has the advantage to be well adapted to quantum graphs since the 
spectra tend to distribute themselves uniformly along the wave number /c-axis, which is not the case along the energy 
axis. The method to obtain the spectra is to consider the rational values of the magnetic flux <i> — p/q, with p,q ^ Z. 
The rational numbers are dense in the real numbers so that we may hope to display the structure of the spectrum by 
plotting the spectra for all the rational numbers <& = p/q up to a large enough integer q. 

If the magnetic flux is rational ^ = p/q, we may assume that the functions satisfy the periodic boundary 

conditions 



("+9^1) (x) = u("^i)(a;) 



(22) 



In this case, the Bloch parameters take their values in the ranges k^ G ^] and ky £ [~f"5 f"] which delimit 

the first Brillouin zone. 



As a consequence of Eqs. (11 1, (12), and (22), the coefficients (13) satisfy the periodic condition 

9n+q gn 



(23) 



and the eigenvalue equation (19) can be solved by requiring that the corresponding characteristic g x g determinant 
is vanishing. 



B. Quantum graphs without magnetic field 



If the magnetic fiux vanishes $ = 0, the eigenvalue equation ( 16 ) can be solved by taking gn+i = gn- The eigenvalues 



are thus given by the zeros of the following function of the wave number k : 

f{k) = (cos klx — cos k^Lj.) sin kly + (cos kly — cos kyly) sin kl^ = (24) 

For fixed values of Bloch's parameters k^ and ky, the spectrum is discrete. A continuous band spectrum is obtained 
by varying Bloch's parameters in the first Brillouin zone delimited hy kx G [—^, f-] and ky G p]- An example 

of band spectrum is depicted in Fig. [l] 

For fixed values of the Bloch parameters, the spectrum is discrete as aforementioned. In this case, we can study the 
statistics of the level spacings : 

[ki+i - ki) 

where ki with i e N are the roots of Eq. (24) : f{ki) = 0. We observe in Fig. [2] that the level spacing distribution 
is empty around zero spacing for kxlx — 1 and kyly — 0.7. This gap is due to the fact that there are only two 
incommensurate lengths and ly in the 2D graph. As a consequence, it is known that the spacing distribution should 
generically present such a gap [6j . 
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Fig. 1: Band spectrum of the 2D graph without magnetic field $ = 0. The bond lengths are l^^l and ly = {l + V5)/2. The 
bands are depicted in the plane of the wave number k — \f2E versus the Bloch parameters k^lx = kyly. 
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Fig. 2: The level spacing probability density P(5') for the 2_D graph without magnetic field "1> = at the values fc^^/a; = f and 
kyly — 0.7 of Bloch parameters. The bond lengths are = \ and ly — (\^ \/5)/2. We observe the expected gap around zero 
spacing. 
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In the classical limit, point particles move with the velocity v on the one-dimensional bonds and are scattered 
stochastically at each vertex [8]. In the case where r = 1, the probability to be scattered in one of the four directions 
of the lattice are equal to 1/4. Accordingly, the particle undergoes a diffusive random walk on the lattice, the properties 

of which can be calculated with the methods of Rcf . [8 . The diffusion coefficients take the values — 2(1+1 ) ^^^^ 

Dy — 2{r~+r^- ^"^^ Kolmogorov- Sinai entropy per unit time is equal to /iks — fTF its positivity shows that the 
classical motion is chaotic. 



C. Quantum graphs with magnetic field 

We solve numerically the generalized Harper equation in order to obtain the energy spectrum as a function of 
the magnetic flux. Several spectra are shown in Figs. |3|a)-(d) in terms of the wave number k — \J2E for the cases 
Ix — 1, 1.5, 2, 2.5 with ly — \ and r —\. 

In the case = ly = ^ and r = 1, the lattice has the square symmetry and the spectrum in Fig. pi a) resembles the 
Hofstadter butterfly [TJI. We can show that it is identical to the Hofstadter butterfly up to a deformation. Indeed, 



Eqs. (17 1 and (18 1 give A = 1 and £ = 4cosfc in this case. Accordingly, Eq. ([19| reduces to Harper's equation which 
leads to the Hofstadter butterfly represented in the plane of the magnetic flux $ versus the energy £. Figure [sf^ a) 
depicts the butterfly versus the wave number k = arccos(£/4). The butterfly is thus only deformed by this change of 
variable. 

However, the anisotropy enhanced by the difference Ix ^ ly closes many gaps [dark zones in Figs. [3|^b)-(d)], which 
reminds a similar phenomenon observed in the anisotropic Hofstadter model [17j . The spectrum of the latter model is 



described by the Harper equation (19) and is darkened when the anisotropy ratio A = |^ differs from unity. We have 
a similar phenomenon in the spectrum of quantum graphs although it is complicated by the common dependence 
of A and £ on the wave number k. We observe that the dark zones appear away from the values ks of the wave 



number for which the anisotropy coefficient is close to unity : A — \ f | ~ 1. The reasons are that the generalized 



Harper equations (19 1 and (20 1 are known to admit extended eigenstates if A is different from unity [2P, and that 
continuous spectra are typically associated with extended eigenstates. In contrast, the spectrum displays fractal 
structures reminiscent of the Hofstadter butterffy around the special values ks where the anisotropy ratio approaches 
unity. 

In the case = 1.5 depicted in Fig. pfb), the anisotropy ratio is unity if |sinfcs| = | sin(3fcs/2)|. This condition is 



satisfled at kg = ... We clearly see in Fig. [3|^b) that the spectrum displays the fractal structures of the 

Hofstadter butterfly around these values, while it darkens away. 

Similarly, we see in the case Ix — '^ depicted in Fig. ISlc) that the spectrum looks locally as the Hofstadter butterfly 
around the special values A:s=f,^,^,^,--- where = 1. 

In the case Ix = 2.5 depicted in Fig. [sj^d), this also happens around the special values kg — 
0.898,1.795,2.094,2.693,3.590,4.189,4.488,5.386,... where A = 1. Here, the situation is more subtle because, for 
instance, the values kg = 1.795 and kg = 2.094 are very close to each other and the anisotropy ratio remains close to 
unity between these values. This explains that the spectrum looks as a Hofstadter butterfly only once in this interval 
and similarly for 4.189 < kg < 4.488. 

We notice that the spectra are periodic in the wave number if the two lengths Ix and ly are commensurate. They 
form non-periodic structures as the wave number increases if the lengths are incommensurate. We flnally note that 
for Ix = ly = n, where n is an integer, the spectrum represents exactly n Hofstadter butterflies for k S [0,7r]. In the 
anisotropic situations depicted in Fig. |3jc) in which Ix — n and ly = 1, we flnd exactly n butterfly-like structures for 
k e [0,7r]. 



III. THREE-DIMENSIONAL QUANTUM GRAPHS 



A. The eigenvalue equation 

The insertion of a third dimension in quantum graphs submitted to magnetic fields is expected to yield drastic 
changes in the spectral properties of these systems. Recent works have been investigating the effects due to three- 
dimensionality in the Hofstadter model [181 |T9j [20] . These studies led to the conclusion that Hofstadter butterffy-like 
fractal patterns still exist in 3D space systems under certain conditions. 

We here show that fractal patterns also persist in 3D quantum graphs. Each bond of the 3D lattice is characterized 
by a vertex coordinate {n,m,l) and by a direction b — 1,2,3, and will be labeled (n,m,l,b). We suppose that the 
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Fig. 3: Spectra of the quantum graphs for the cases : (a) Ix = 1, (b) Ix = 1-5, (c) Ix = 2, and (d) Ix = 2.5 with ly — 1 and 
r = 1. The wave number k = y/lE is along the horizontal ajcis and the magnetic flux $ along the vertical axis. 



magnetic field is still oriented along the z-axis and that we work with the Lorentz gauge where A — (0,-Ba;,0), so 
that the wave function along the bond (n, m, 3) is a solution of the Schrodinger equation 



(26) 



We can then treat this system using the method described in the previous Sec. [TT] We impose boundary conditions 
similar to the ID case with 

^^(n,m,l,l) ^^(n-l,rn,l,l) ^^(n,m,l,2) ^^(n,m-l,i,2) 

■(^ = 0) iz, + — i:. — (y = 0) 1:. ^ ^y) 



dx 
i B nl„ 



dx dy 



dy 



+ 7 



= 0) — {z = 



dz 







where the parameter 7 S [0, 1] maps the 2D system (7 = 0) to the 31? isotropic system (7 = 1). 
The wave function is explicitly written as 



sm Kt^ 



sin kz 



n.m,/ + l . , , 

smfct^ 



Setting 

with kx € ky e [^f-i f-] and fc^ G [-~f-,f-], one finds the generalized 3D Harper equation 



(27) 

(28) 
(29) 



2 A cos(27r<i>7i - kyly) + 2 A cos kj 



gn ^ £gn 



with the coefficients 



A = 



sin klr, 



sin kL, 



A = 7' 



2 sm k^x 



f = 2 (cotg fc^a; + cotg /cZj^ + 7^cotg /cZ^) sin fc/a; 



(30) 

(31) 

(32) 
(33) 
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Fig. 4: Band spectrum of the 3D graph without magnetic field $ = 0. The bond lengths are = l,ly^{l + V5)/2, and 
Iz = exp(f). The bands are depicted in the plane of the wave number k — \f2E versus the Bloch parameters k^lx = kyly — kzh- 

As in the 2D system, these coefficients all depend on the wave number k. 

B. Quantum graphs without magnetic field 

If the magnetic flux vanishes $ = 0, we may assume that the solution is periodic with (?„+i = Qn, in which case the 
eigenvalue equation (30 1 with the coefficients (31)-(33l becomes 



f{k) — (cos kLjr — cos fca-Zj.) sin kly sin kl^ + (cos kly — cos kyly) sin kl^ sin kl^ + 7^ (cos kl^ — cos k^lz) sin klr^ sin kly — 

(34) 

Because of its spatial periodicity, the system has continuous band spectra, an example of which is depicted in Fig. 
|4] Contrary to the 2D system, Wigner repulsion manifests itself in the level spacing statistics. The Wigner repulsion 
has for consequence that a typical level spacing probability density behaves as P{S) ^ S" for S" ^ 0, as observed in 
Fig.jSj Consequently, there is no gap at small level spacings contrary to what happens in Fig. [2] for the 2D case. The 
reason is that the 3Z) lattices typically have three incommensurate lengths, which is the minimum number for Wigner 
repulsion to manifest itself [6]. 



C. Quantum graphs with magnetic field 



We solved the generalized ZD Harper equation (30) with <f> 7^ in order to find the spectrum associated to the 
system for different values of 7. As shown in Fig.^6f the butterfly seems to lose its initial symmetric shape as 7 
differs from zero, and eventually forms a new fractal structure for 7 = 1. It is worth noticing that the many gaps 
which compose the Hofstadter butterfly are conserved as one maps the 2D system to the 3D system. Figure |6] depicts 
the spectra corresponding to a single value of the Bloch parameter kz, which corresponds to the situation of a very 
flattened 3-D system. As the graph thickens, kz takes more values between —i^jlz and tt/^z and the spectrum appears 
as the superposition of several deformed butterflies. As the three-dimensionality becomes important, gaps close and 
the spectrum darkens as shown in Fig. [Tj 
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Fig. 5: The level spacing probability density P{S) for the 3D graph without magnetic field $ = at the values k^lx ~ 1-5, 
kyly = f and k^l^ = 0.7 of Bloch parameters. The bond lengths are Ix = 1, ly — {1 + \/5)/2, and Iz = exp(l). We observe the 
behavior characteristic of Wigner repulsion : P{S) ^ S for 5^0. 




k k 

Fig. 6: Spectra of the quantum graphs for : (a) 7 = 0.5, (b) 7 = 1, with Ix ~ ly = Iz = 1 and kz = 0. The wave number 
k — \/2-E is along the horizontal axis and the magnetic flux $ along the vertical axis. 



IV. QUANTUM HALL EFFECT ON QUANTUM GRAPHS 



A. Kubo formula on quantum graphs 



In order to study the integer quantum Hall effect, we consider independent Fermions on a 2D quantum graph 
and evaluate the antisymmetric component (Jxy — —(Jyx of the conductivity tensor with Kubo's formula for zero 
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Fig. 7: Spectra of the quantum graphs for 7 = 1, with Ix = ly = Iz ~ 1 and = nyr, n G Z. The wave number k — \/lE is 
along the horizontal axis and the magnetic flux $ along the vertical axis. 

temperature. The latter relates Gxy to the current intensity J — {Jx, Jy) and can be written in the following way 



iV 



N 



(35) 



where |N) denotes the eigenstates of the many-particle Hamiltonian H of eigenvalue E-^, |0) is the ground state, and 
V = J2(n m b) is the system "volume" given by the total lengths of the bonds in a unit cell of the periodic lattice. 
lb denotes the length of the bond b. 

We define the current density which circulates along the bond (n, to, b) in the direction x by 



1 / J^(n,m,6) ja,(n,m,b)t \ 
j{n,m,b) _ I ||^(n,m.b)t ^(n,m,z) j ^(n,m,6)^(n,m,6)t^(n,m,fc) 

^ 2i \ dx dx I ^ 



(36) 



where A^"'™'''^ is the x-component of the potentiel vector along the bond (n, m, b) and ^(■^'"^■'>) is the field operator 
defined on the bond. In the second quantization formalism, the field operator and the adjoint field operator are 
respectively given by 



(n,m,6)* -.-j- 



(37) 
(38) 



where and aj, are the annihilation and creation operators satisfying the anticommutation relation : ci^a^^, +a\,di, = 
6vu' ■ If the system is in the Fock state |N) — \NiN2N^ ■ ■ ■ Ni, ■ ■ ■) with N^, = or 1 Fermion on each single-particle 
wave function ijju, the total energy of the system takes the value 



(39) 



where ev is the energy corresponding to the single-particle wave function -01/ and is the corresponding occupation 
number {N^, = 0, 1). 

The current intensity which circulates in a unit cell of the lattice is defined by 



(n,m,fc) 



(40) 



B Quantization of transverse conductivity 
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where the sum over (n, m, b) extends over all the bonds of the unit cell. Substituting the expression (36) of the current 
density and using Eqs. (|37[)-(38 1, we find that the current intensity is given by 



(n,m,b) v.v' 



1 / J I (n,7n.b) 1 I {n.7n,b)* \ 

n,m,b) ^(n,m,b)* ^(n,m,b) 



2i 



dx 



dx 



K. K (41) 



(42) 



The quantities (y'\vx\v) which are here introduced are the matrix elements of the particle velocity in the single-particle 
wave functions V'y • The scalar product in the space of the single-particle wave functions on the graph is defined by 



/ i.\ Jo 
{n,rn,b) 



(43) 



where the sum J2(nmb) accounts for the contribution of all the bonds in the unit cell. 

Since the operator contains one annihilation and one creation operator, the matrix element (N| j^|0) is non 
vanishing only for Fock states 



|N) = I111---111 HI- ••111 I 000- ••000 1 OOO^^^) 



(44) 



with one hole and one particle, ep denotes the Fermi energy. Indeed, applying an annihilation operator followed by a 
creation operator on the ground state gives such a state up to a phase : 



4aa|lll---lll I OOO^^^OOO^^^) = (-1)°+'^^|111^^^111 lll^^^lll I OOO^^^OOO 1 000 • 



(45) 



where N-p is the total number of Fermions and a the integer labelling the corresponding single-particle state. Therefore, 
the matrix element of the current intensity is given by 

(N|J.|0) = (-l)"+^-(/3|fi.|a) (46) 

for the Fock state (44l and zero otherwise. Moreover, we notice that the energy of this Fock state (44l is equal to 

Ej^-Eo^ep- e„ (47) 

The energy of the single-particle state |a) is below the Fermi energy < ep, while the situation is opposite for the 
state 1/3) : e/3 > ep. 

Therefore, the conductivity (35 1 can be expressed in terms of single-particle wave functions of the quantum graph 

1 



as 



iV 



E E 

eQ<eF «/3>eF 



'\vx\a){a\vv\P) -{x ^y) 

(ea - 6/3)2 



(48) 



B. Quantization of transverse conductivity 



We consider the single-particle wave functions (111 and (12 1 with the periodic boundary conditions (22). The 



corresponding single-particle Hamiltonian is defined by 



1 2 1 2 

-H" = 2 - Ax + K) + - {-idy - Ay + ky) 



(49) 



where and ky are the components of the wave vector. The latter are Bloch parameters which take their values on 
the first Brillouin's zone of the reciprocal space : k^ — r^, ^] and ky E [— f^, t^]. The velocity operator is given in 
terms of this Hamiltonian by 



dH 
dkr 



(50) 
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We consider the Hamiltonian eigenstates Uq"''''' which satisfy Schrodinger's equation 



Hut'^ = e^ut"^ (51) 



and the periodic boundary conditions ( 22 1 . Differentiating the eigenvalue equation ( 51 1 with respect to one component 
of the wave vector and taking the scalar product with another eigenstate, we get 

BM rirv 

(/3|^.|«>(= (/3|g^l"> - - (52) 

with the short notations |a) = \ua) and \dk^a) = \dk^Ua). 

Using the relation {I3\dk^a) = —{dk^(3\oi) obtained by differentiating the orthonorniality condition {a\(3) = 5ap 
together with the completeness relation (X]e„<eF + Se„>eF) ~ ^'^'^ conductivity can be written 

as [10] 

iV ^^^^\dky^dK/ \dkjdky/ ^ ' 



If the Fermi level falls inside a spectral gap, the sum over all the states \a) below the Fermi level can 

be decomposed into a sum over the occupied bands and a sum over the states \a) inside a band. This latter can be 
performed as an integral over the values of the Bloch parameters in the first Brillouin zone which forms a torus T^. 
Accordingly, we have that 

ea<eF occ. bands in occ.band occ. bands 

Finally, using Stokes theorem in order to transform the two-dimensional integral over the torus into a line integral 
over the border of the torus, the transverse conductivity becomes 

-^. = 2^ E /^"iI^^^'^^ + ^^Im^)*^ ^^^^ 

occ. bands y / 

where the dimensionnal factor /h with h = 2Trh has been reintroduced |10j . 

The above expression has a well-known topological interpretation : it is the integration of Berry's curvature J- = 
(^1 ^)dkf' A dk" over the base space of a principal fiber bundle F{T^,U{1)). The number Chi = 4n Ir^ ^ 
is a topological invariant referred as Chern number and is necessarily an integer. Thouless et al. showed that the 
Hall conductivity (55) is given by a^y — J2e^<eF i^r — tr-i) for Hofstadter's model, where tr is the solution of a 
diophantine equation 

r ^ q Sr + ptr (56) 

which gives the r*'^ gap's position of the Hofstadter spectrum in terms of the integers p and q of the magnetic flux 
<i> — p/q. When the Fermi level is exactly situated in the r"' gap, one finds that a^y — t^e'^/h. The Chern number 
being an integer one gets the quantization law observed by von Klitzing [5]. 

Now, we show that this result extends to quantum graphs. We suppose that the wa ve functions accumulate the 
phase Ua ~ exp(z6'a) along the border OT^ of the torus where the line integral of Eq. (55) is carried out, in which case 
one can write ( 55 ) in the following way 



=2 



X- I fdO^^k^^^dky] (57) 



y 



In the weak-coupling limit r ^ 0, the anisotropy parameter (17) is large and the generalized Harper equation (19) 
reduces to the approximate eigenvalue equation 

cos kly ~ cos (27r<i>n — kyly) (58) 

with n g Z. These equations form crossing curves in the plane of the wave number k versus the Bloch parameter ky. 
The crossings are exact if r = 0, but avoided crossings exist if r is small but non vanishing. In this case, it is possible 
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to calculate the phase Oa accumulated by the eigenstate Ua at each avoided crossing, in a way similar to the one 
shown for tight-binding models by Kohmoto [16j . We obtain that the eigenstates Uk^k undergo the transformation 



Jqk^l^{tr-tr--i) 



(59) 



when a loop is performed around the first Brillouin zone. One eventually finds that the Hall conductance is given by 



2nh 



E 



^ de 

dkx TT, 



27r(tr-tr-l) 



{tr — tr-l) 



e2 



(60) 



where we have supposed the Fermi level ep in the r gap 



C. Phase diagrams of transverse conductivity 



The phase diagram which describes the integer quantum Hall effect for Hofstadter's model has been introduced by 
Osadchy and Avron (llj . This diagram shows in a beautiful fractal structure the value of the Hall conductivity as a 
function of the magnetic flux and the Fermi energy of the system. Such a phase diagram can be drawn for quantum 
graphs as well. The result is straightforward for the self-dual case Ix = ly [17l|2T], for which the diagram is a nonlinear 
deformation of the one obtained by Osadchy and Avron [11| . In order to get the phase diagram of an arbitrary graph, 
we have solved numerically the generalized Harper equation and resolved the diophantine equation ( 56 ) for each gap 
of the spectrum. The phase diagram for the case = 3 and ly = 1 is drawn in Fig. |8]as a function of the Fermi wave 
number fcp — •\/2eF and the magnetic flux <I>. The quantum phases correspond to the different integer values of the 
Hall conductivity computed inside the gaps. 

A phase diagram describing the IQHE is also obtained for the three-dimensional system described in Sec. The 
Chern numbers evaluated inside the numerous gaps of the initial butterfly (7 = 0) are maintained while the spectrum 
undergoes the transformation 7 — > 1 . We draw the phase diagram corresponding to the 31? quantum graph in Fig. [9] 



V. CONCLUSION 



In this article an important aspect of the quantum transport on 2D and 3-D graphs have been studied namely the 
quantization of the system's Hall conductivity. 

First, we have obtained the energy spectra of quantum graphs without and with magnetic field. We showed that 
their eigenvalue equation can be mapped onto a generalized Harper equation in the case of a 2D rectangular lattice. 
A 3D rectangular lattice has also been considered. 

In zero magnetic field, the graphs have continuous band spectra because of their spatial periodicity and the spectra 
are discrete at fixed values of Bloch's parameters. The 2D and 'iD graphs are shown to differ by their level spacing 
statistics. Indeed, the 2D graph has at most two incommensurate bond lengths so that its level spacing distribution 
typically presents a gap around zero spacing. In contrast, the 31? graph has at most three incommensurate bond 
lengths which is sufficient for Wigner repulsion to manifest itself in the level spacing statistics. On the other hand, 
both the 2D and 31? graphs are classically chaotic with a positive Kolmogorov- Sinai entropy per unit time in the 
classical limit. 

In non-zero magnetic field, we have obtained fractal energy spectra. A deformed version of Hofstadter's butterfiy is 
recovered in the case of a 2D lattice with the C4 square symmetry. If the lattice becomes anisotropic, some gaps are 
filled and the corresponding zones darken in the spectrum due to the appearance of continuous parts in the spectrum. 
Nevertheless, other gaps remain which are characterized by Chern's topological quantum numbers. We show that the 
transverse conductivity is quantized in terms of Chern's numbers, as in the integer quantum Hall effect. We construct 
the fractal quantum phase diagrams of the transverse conductivity. 
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Fig. 8: Spectrum and phase diagram for a 2D quantum graph characterized by Ix = 3, ly = 1, and r — 1. The horizontal axis 
is the Fermi wave number fcp — \f2tY and the vertical axis the magnetic flux $. The different phases correspond to the integer 
values of the Hall conductivity computed inside the gaps. As in the situation depicted in Fig. [sjc) we have ~ n and ly = 1, 
with n = 3, and we find exactly n — 3 butterfly-like structures for k £ [0,it]. 

In conclusion, quantum graphs show rich structures such as fractal spectra and reveal interesting quantum 
properties such as those emphasized in this work. These versatile models are promising for the exploration of 
quantum phenomena in condensed matter systems such as the quantum Hall effects. 
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Fig. 9: Spectrum and phase diagram for a 3-D quantum graph characterized by Ix ~,ly — Iz = 1. The horizontal axis is the 
Fermi wave number fep = y/2eF and the vertical axis the magnetic flux 3>. The different phases correspond to the integer values 
of the Hall conductivity computed inside the gaps. 
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